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ABSTRACT 

The nearby low-luminosity active galactic nucleus (LLAGN) NGC 4258 has a weak radio continuum 
emission at the galactic center. Quasi-simultaneous multi-frequency observations using the Very Large 
Array (VLA) from 5 GHz (A6cm) to 22 GHz (Al.Scm) showed inverted spectra (a > 0, S*^ oc i/") in 
all epochs, which were intra-month variable, as well as complicated spectral features that cannot be 
represented by a simple power law, indicating multiple blobs in nuclear jets. Using the Nobeyama 
Millimeter Array (NMA), we discovered a large amplitude variable emission at 100 GHz (A3mm), 
which had higher flux densities at most epochs than those of the VLA observations. A James Clerk 
Maxwell Telescope (JCMT) observation at 347 GHz (A850 /im) served an upper limit of dust con- 
tamination. The inverted radio spectrum of the nucleus NGC 4258 is suggestive of an analogy to our 
Galactic center Sgr A*, but with three orders of magnitude higher radio luminosity. In addition to 
the LLAGN M 81, we discuss the nucleus of NGC 4258 as another up-scaled version of Sgr A*. 

Subject headings: galaxies: active — galaxies: individual(NGC 4258) — galaxies:jets — galaxies: 
Seyfert — radio continuum: galaxies — submillimeter 



1. INTRODUCTION 

The nucleus of NGC 4258 (7.2 Mpc, Herrnstein et 
al. 1999) is one of the three closest radio continuum 
sources as extremely sub-Eddington super-massive black 
hole (SMBH) systems, following Sgr A* (8 kpc, Reid 
1993; Eisenhauer et al. 2003) at the Galactic center 
and the nucleus of M 81 (3.6 Mpc, Freedman et al. 
2001), as shown in Table 1. Their low- luminosities are 
believed to be caused by the radiative inefficiency of 
optically-thin accretion inflows in low accretion rates. An 
advection-dominated accretion flow model (ADAF) can 
fairly successfully explain the observed broadband spec- 
tra of Sgr A* (e.g., Narayan et al. 1995; Manmoto et al. 
1997). In this model, synchrotron emission of thermal 
electrons is responsible for a highly inverted radio spec- 
trum. However, this model has two problems: (1) The 
observed slightly inverted spectrum means the presence 
of a significant radio excess compared to the ADAF spec- 
trum. The radio excess requires an additional emitting 
component. A model of jet -I- ADAF has been proposed 
to explain the radio excess (Falcke & Biermann 1999; 
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Yuan et al. 2002). In another model, an inflow of hybrid 
plasma consisting of both thermal and nonthermal par- 
ticles is adopted for Sgr A* (Ozel et al. 2000). (2) The 
upper limit of the integrated electron density of accreting 
plasma, constrained by the detections of radio polariza- 
tion in Sgr A* (Bower et al. 2003), is much lower than 
that predicted by the ADAF model. A sufficient mass 
loss is required before advecting across the event horizon. 
A radiatively inefficient accretion flow (RIAF) model was 
also proposed, involving a substantial mass-loss outflow 
and a non-thermal component of electrons (e.g.. Yuan et 
al. 2003). Thus, an outflow or a jet may be essential for 
modeling the radio spectrum. However, the presence of 
an outflow or a jet in Sgr A* has never been established 
in any observational investigation. 

Sgr A* shows a mysterious inverted spectrum (a > 0, 
Si, oc i/+") over the radio regime (e.g., Falcke et al. 
1998; An et al. 2005), which cannot be explained by the 
simple ADAF model. The intrinsic size is ~ 1 AU at 
86 GHz, corresponding to ~ 13 times the Schwaltzschild 
radius (Shen et al. 2005); there has been no evidence for 
any radio structure of jet or disk. The M 81 nucleus 
also shows an inverted spectrum similar to Sgr A* but 
is a much more powerful radio source; therefore, a so- 
called "a scaled-up version of Sgr A*" (Renter & Lesch 
1996), or "a bridge source" to help span the gap be- 
tween Sgr A* and other low-luminosity active galactic 
nuclei (LLAGNs) (Markoff et al. 2008), has been sug- 
gested. The M 81 radio source shows a slightly resolved 
jet-like structure, with a size of 700x300 AU at 22 GHz 
(Bietenholz et al. 1996). However, the jet-like compo- 
nent in M 81 makes only a minor contribution to the 
total flux; the substantial contribution to the observed 
inverted spectrum is from an unresolved core. 

NGC 4258 is a LLAGN of type 1.9 Seyfert galaxy (Ho 
et al. 1997) with a weak radio continuum component at 
its galactic center. Total flux densities at 22 GHz and 
1.5 GHz were ^ 3 mJy and 0.5 mJy, respectively, in very- 



2 



Doi et al. 



TABLE 1 
The nearest three active nuclei. 



Name D Ref. Mbh Ref. P22 &22 Rcf. Lx Lx/Ljsdd Rcf. 

(Mpc) (Mq) (W Hz-1) (mas) (AU) (Rs) (ergs s-i) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 

NGC 4258 7.2 a 3.9 x lO'' a 2 x lO^^ 0.4 3000 3700 g 8.0 x 10*0 1.6 x 10-^ j 

M 81 3.6 b 7.0 x lO'' e 1.9 x lO^" 0.18 650 460 h 1.5 x 10"° 1.7 x 10"^ k 

Sgr A* 0.008 cd 4 x 10^ f 5.7 x 10^^ 0.46 4 50 i 2.4 x 10^3 5 x lO'^^ 1 

Note. — Col. (1) source name; Col. (2) distance; Col. (3) reference of distance; Col. (4) black hole mass; Col. (5) reference of black 
hole mass; Col. (6) radio power at 22 GHz; Col. (7) projected size of radio source at 22 GHz in milli-arcsecond. For NGC 4258, the offset 
of a northern nuclear jet from a dynamical center. For M 81, the major axis of the fitted size of the nuclear component. For Sgr A*, 
an intrinsic source size estimated by calibrating interstellaj: scattering, see references in detail; Col. (8) size of radio source at 22 GHz in 
AU; Col. (9) size of radio source at 22 GHz in Schwarzschild radius; Col. (10) reference of radio size; Col. (11) absorption-corrected X-ray 
luminosity at 2—10 keV; Col. (12) ratio of X-ray to Eddington luminosities. Col. (13) reference of X-ray luminosity. 

References. — a; Herrnstein et al. 1999, b: Freedman et al. 2001, c: Reid 1993, d: Eisenhauer et al. 2003, e: Devereux et al. 2003, f; 
Gillessen et al. 2009, g: Herrnstein et al. 1997, h: Bietenholz et al. 1996, i: Shen et al. 2005, j: Yamada et al. 2009, see also Pruscione et 
al. 2005, k: Markoff et al. 2008, see also Ishisaki et al. 1996, 1: Baganoff et al. 2003. 



long-baseline interferometry (VLBI) observations (Her- 
rnstein et al. 1997; Cecil et al. 2000). Multi-frequency 
observations at 0.5" resolution toward the nucleus re- 
sulted in 2.1, 1.9, and 2.8 mJy at 5, 8.4, and 15 GHz, 
respectively (Nagar et al. 2001). Thus, it remains possi- 
ble that the NGC 4258 nucleus has an inverted spectrum 
similar to Sgr A* and M 81 nucleus. 

The VLBI imaging showed an apparently small jet with 
a projected scale of ^ 2 x 10^ AU at 22 GHz, upwcUing 
from the dynamical center of a Keplcrian disk rotating 
around a dark mass of 3.9 x 10^ Mq and well established 
by water maser observations (Miyoshi et al. 1995; Her- 
rnstein et al. 1999). The putative accretion disk on the 
SMBH has an almost cdgc-on orientation with the rota- 
tion axis pointing almost north {PA = —8°) and tipped 
82° to the line-of-sight. The jet closely aligns in projec- 
tion on the sky with the rotation axis. The jet extends 
up to ~ 5 X 10^ AU at 1.5 GHz in a VLBI image (Ce- 
cil et al. 2000), which may be physically connected with 
the kpc-scale jet-like structures ("anomalous arms": e.g., 
van dor Kruit et al. 1972; Hyman et al. 2001). Thus, the 
key difference between NGC 4258 and Sgr A* and M 81 
is that the jet in NGC 4258 has definitely been identi- 
fied as the radio emitting origin by VLBI imaging. The 
peaks of north-south jets are significantly shifted from 
the putative position of SMBH, implying that the ma- 
jor fraction of total flux is from jets rather than from an 
innermost accretion flow. If the nucleus of NGC 4258, 
just like that of M 81, is an another scaled-up version of 
Sgr A*, the radio emitting processes of the mysterious 
inverted-spectrum radio emission of Sgr A* and M 81 
could be responsible for jets. However, an inverted spec- 
trum over the radio regime has not yet been established 
for the NGC 4258 nucleus. 

In the present paper, we report the continuum spectra 
and flux variability of radio emission of the NGC 4258 
nucleus at centimeter to submillimeter bands. Observa- 
tions and data reductions are described in § 2. Results 
are presented in § 3, and their implications arc discussed 
in § 4. Finally, we summarize in § 5. At a distance of 
7.2 Mpc to NGC 4258, 1" corresponds to 35 pc. 

2. OBSERVATIONS AND DATA REDUCTIONS 

We reduced raw data from the archival data of the 
James Clerk Maxwell Telescope (JCMT) on Mauna Kea 
at 347 GHz, the data of our observations conducted 



with the Nobeyama MilHmeter Array (NMA) at the 

Nobcyama Radio Observatory (NRO) at ^ 100 GHz, 
and from the archives of the Very Large Array (VLA) 
at the National Radio Astronomy Observatory (NRAO) 
from 4.9 to 43 GHz. 

2.1. JCMT data 

We retrieved archival data for observations of 
NGC 4258 on March 18, 1998 using the SubmilHmetre 
Common User Bolometer Array (SCUBA) (Holland et 
al. 1999) on the JCMT at 347 GHz (850 ^um) in jiggle- 
map mode. The SCUBA covers a field of view of ~ 2', 
which is two times larger than that covered by the NMA. 
The data were processed with the SCUBA User Reduc- 
tion Facility (SURF) package. We applied the standard 
reduction procedures including flat fielding, flagging of 
transient spikes, correcting the extinction, pointing cor- 
rection, sky removal and flux density scaling. An Uranus 
scan gave us a flux scaling factor, with an uncertainty as- 
sumed to be ^15%, and an effective beam size of 15.1", 
corresponding to ~530 pc. The measured R.M.S. of the 
image noise was 10.5 mJy beam~^ on blank sky. 

2.2. NMA observations 

We observed the nuclear regions at ^ 100 GHz 
(A3 mm) using the NMA. Array configurations were AB, 
C, and D with maximum baseline lengths of 351, 163, 
and 82 m, respectively. Observations were originally con- 
ducted for HCN( J = 1-0) and HCO+( J 1-0) molecu- 
lar emission lines, except for observations conducted only 
for continuum between May 13 and 15, 2005. Visibility 
data were obtained with double-sided receiving systems. 
In the data for molecular-line observations, we use only 
line-free data from the upper-side band at a center fre- 
quency of 100.777 GHz with a bandwidth of 1 GHz (the 
emission lines were in the lower-side band at a center fre- 
quency of 88.777 GHz). In the data for continuum obser- 
vations, we use two sidebands at 89.729 and 101.729 GHz, 
both of which were line-free. The typical system noise 
temperature was about 140 K in the double-sided band. 
For gain calibration, we made scans between the target 
and a reference calibrator, 1150-1-497, every 20 or 25 min. 
Flux scales of the calibrator were derived by relative com- 
parisons with the known flux calibrators such as Uranus 
that were observed quasi-simultaneously at almost the 
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same elevations in order to avoid a residual error of cal- 
ibration from differential atmospheric attenuation. The 
uncertainty in the absolute flux scale is expected to be 
less than 15%. A scan of around 30 min was performed 
for bandpass calibration using bright quasars. 

The data were reduced using the UVPRDC-II pack- 
age (Tsutsumi et al. 1997) by standard manners, includ- 
ing flagging bad data, baseline correction, opacity cor- 
rection, bandpass calibration, and gain calibration. For 
the data observed in May 2005, we combined the upper 
and lower sidebands to obtain a higher sensitivity; the 
resultant center frequency was 95.729 GHz. Each daily 
image was synthesized from visibilities in natural weight- 
ing and deconvolved using the IMAGE task of the Astro- 
nomical Image Processing System (MPS) developed at 
the NRAO. FWHMs of the synthesized beams were typ- 
ically ^ 3", 5", and 7" in the AB-, C-, and D-array con- 
figurations, respectively. These resolutions correspond to 
-100, 170, and 240 pc. The R.M.S. of image noise was 
estimated from statistics on blank sky using the IMEAN 
task of AIPS and was used to determine 3cr upper limits 
in the cases of negative detection. 

2.3. VLA data 

We retrieved 12 archival data sets from VLA. Our se- 
lection criteria for the data were (1) observation toward 
the nucleus in continuum mode, (2) multi-frequency ob- 
servation quasi-simultaneously at three or four bands in a 
range of 4.9-22 GHz or observation including at 43 GHz, 
and (3) observation with high-spatial resolution (< 1.2") 
of A- or B-configurations to extract only a nuclear com- 
ponent. We reduced the data using the AIPS by following 
a standard procedure. 3C 286 scans served as flux scaling 
factors, with an uncertainty of 3 and 5% that we assumed 
at <22 GHz and >22 GHz, respectively. 

In data obtained on January 04, 1997, IF2 at 15 GHz 
was very noisy, so we used only data of IFl. In data 
obtained on March 06, 1997 and September 05, 1998, 
22-GHz observations had been conducted to set IPs to 
three frequencies, corresponding to the well-known mega- 
masers of the low-velocity component near the systemic 
velocity and of high- velocity components red-shifted and 
blue-shifted relative to the systemic velocity, respectively. 
However, these observations were done in continuum 
mode. We used only the blue-shifted maser's band, at 
which the maser emission is sufficiently weak (Nakai et 
al. 1993) to not affect our continuum photometry. 

All target images were synthesized with natural 
weighting and deconvolved with the IMAGR task of AIPS. 
The resultant FWHM of the synthesized beam is 0.1- 
2". The typical R.M.S. of image noise estimated from 
statistics on blank sky using the IMEAN task of AIPS was 
0.45, 0.10-0.20, and 0.05-0.10 at 43 GHz, 22-15 GHz, 
and 8.4-4.9 GHz, respectively. 

3. RESULTS 

The results of the flux measurements are listed in Ta- 
ble 2, and the details are described in the following sub- 
sections. An average of flux densities over all epochs at 
each frequency, a reduced chi-square, and a probability 
of no variability are also shown at the lowest three lines 
in the table. Flux variability is found at all frequencies 
at greater than 95% confidence, except for 43 GHz and 
347 GHz, where the number of epochs was inadequate 
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Fig. 1. — All measurements for the nucleus of NGC 4258 are 
shown in Table 2. Upper limits at ~100 GHz are represented by 
downward arrows. The 347-GHz measurement is the intensity of 
the nuclear region in the JCMT beam, which can be interpreted 
as the upper limit of nuclear emission. Dashed line represents the 
upper limit of dust continuum of gray body emission with /9 = 1 
(i.e., a = 3). The Scr upper limit of flux density of continuum 
emission at 225 GHz observed on May 4, 2004, which was measured 
with a beam size of 3.0" X 2.0" (Sawada-Satoh et al. 2007), smaller 
than the beam size of the NMA, is shown as a reference. 

for statistics. At 100 GHz, six upper limits are included. 
Statistics for the NMA flux densities are performed using 
only detected data. The spectral plot of all the data is 
shown in Fig. 1, which apparently shows a global trend 
of an inverted spectrum over the radio bands. 

The JCMT image is shown in Fig. 2, which is the only 
case in which some radio structure is seen in our images. 
The light curve plot for measured NMA flux densities is 
shown in Fig. 3. The spectral evolution obtained by the 
VLA observations is shown in Fig. 4. 

3.L The upper limit of dust contamination 

The JCMT image at 347 GHz (A850 /xm) of the nuclear 
region (~ 4.3 x 4.3 kpc; Fig. 2) shows a structure elon- 
gated along the major axis of the host galaxy and with 
no nuclear concentration. Such a structure is very simi- 
lar to that of CO emission (Regan et al. 2001; Heifer et 
al. 2003; Sawada-Satoh et al. 2007). These suggest that 
most of the submillimeter emission originated not in the 
AGN continuum but in interstellar dust associated with 
the host galaxy. The flux density from the active nucleus, 
at the point indicated by a cross in Fig. 2, contributes 
less than the peak intensity. We averaged intensity over 
the JCMT beam of 15.1" centered at the nucleus in or- 
der to obtain an upper limit of the flux density of AGN 
emission. The upper limit may be contaminated by the 
CO (3-2) emission of 1 mJy in the SCUBA bandwidth 
of ~ 40 GHz, considering the velocity-integrated inten- 
sity peak of the CO (1-0) map (Heifer et al. 2003) and 
the intensity ratio CO(3-2)/CO(l-0) (Mao et al. 2010). 

In the present study, the JCMT measurement does 
nothing more than providing the upper limit of AGN 
continuum emission contaminated by dust emission, as 
shown by a dashed line in Fig. 1. The gray body spec- 
trum of dust emission is S'^^^^ oc in the Rayleigh- 
Jeans regime at less than about 1 THz, where /3 is the 
dust emissivity, which is believed to be between 1 and 2 
(Hildebrand 1983); the fi value has been determined to 
be 1.3 ± 0.2 by JCMT observations for nearby galaxies 
(Dunne et al. 2000). In this paper, we adopt /? = 1, which 
is the most influential case in the spectral extrapolation 
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TABLE 2 

Flux densities of radio continuum emission at the nucleus of NGC 4258 and spectral indices. 



uate 


Array 


■is 


■58.4 


5'l5 


522 


•343 


Q 

•J 100 


•3347 


^22 
•^5 




config. 


(mjy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 




w 








(5) 


(6) 


v) 






(10} 


U7 Dec lyyo 


VIjA-A 


l.SitO.l 


1.9it0.1 


2.2±0.2 


2.4±0.3 








U.21±U.U4 


22 Dec 199d 


VLA-A 


2.0±0.1 


2.2it0.1 


2.2±0.2 


2.2±0.2 








0.0d±0.0d 


29 Uec lyyb 


VLA-A 


1.9±0.1 


2.3±0.1 


2.7±0.2 


2.5±0.3 








0.21±0.07 


61 uec 199d 


VLA-A 


1.9±0.1 


2.3±0.1 


2.7±0.3 3.1±0.4 








r\ oo 1 r\ OI 

U.cJ2±U.Ul 


U4 Jan 1997 


TTT A A 
V1jA-A 


1.8it0.2 


2.0it0.1 


3.5±0.2 










U.d4±U. 2a 


Uo Jan lyy i 


■\ /"T A A 
V LA- A 


2.3it0.1 


2.3it0.1 


2.6±0.3 


3.6±0.4 








U.19±U.12 


Ub Mar 1997 


VLA-B 


S.litO.S 


2.7±0.2 




5.4±0.3 








r\ £;o 1 r\ 1 o 

0.b2±0.19 


19 Mar 1997 


VLA-B 


2.7±0.3 


2.8±0.3 




3.2±0.4 








r\ 1 r\ 1 r\ 01 

U.1U±U.U1 


zo rCD 199(5 


■\ /"T A A 










7.1it0.9 








lo Mai lyyo 


JOM 1 














9o. / iL2U.f5 




Uo bcp 199o 


V LA-ri 


2.5it0.3 


2.2it0. 1 


3.6±0.4 3.0±0.4 








n oi^_Ln OI 
U.2o±U.21 


Uo bcp 1999 


V LA- A 


2.1it0.2 


l.SitO.l 


2.5±0.3 










r\ 1 fi_Lr\ O/t 
U.19±U.24 


zi Uec zUUU 


■\7'T A A 
V1jA-A 




Z.DltU.Z 














21 Uec zUUl 


INMA-U 












10. Oil. 7 






21 Mar zUU2 


iNMA-O 












18.1di4.0 






22 Uec 2UUo 


■\TTV/r A 
INMA-U 












5.2dil.2 






14 Jan 2UU4 


INMA-ArS 












< 5.9 






31 Mar 2004 


NMA-D 












6.2ib2.0 






01 Apr 2004 


NMA-D 












< 8.3 






Uo Apr 2UUo 


IN MA-U 












4.6dil.8 






U8 Apr 2UUo 


TV T A :T A 

INMA-C 












8.0ib3.3 






±o ivioiy zuuo 


IN W\.r\- lJ 


















14 May 2005 


NMA-D 












10.4±4.0 






15 May 2005 


NMA-D 












7.0±2.1 






25 Mar 2006 


NMA-D 












< 9.4 






29 Mar 2006 


NMA-D 












< 9.4 






31 Mar 2006 


NMA-D 












< 7.1 






01 Apr 2006 


NMA-D 












< 8.6 






Average 




2.0±0.1 


2.2±0.1 


2.7±0.2 


3.3±0.5 


6.9±0.2 


7.1±0.9"^ 




0.21ib0.04 






4.11 


5.66 


4.73 


30.33 


0.06 


1.97"^ 




25.81 






< 10-4 


< 10-4 


< io-» 


< 10"" 


0.804 


0.046"= 




< 10-4 



Note. — Col. (1) observation date; Col. (2) telescope, and its array configuration; Col. (3)-(9) fiux densities; Col. (10) spectral index, 
derived by least-squares fit to Su cx at centimeter wavelengths. 

^Intensity of the nuclear region in the JCMT beam as an upper 
limit of nuclear emission. 

''Nagar et al. (2001) 

"^Derived from only nine detected data sets. 

to lovi^cr frequencies. As a result, the dashed hue in Fig. 1 
is useful to estimate the upper limit of dust contamina- 
tion in the millimeter photometry of the nucleus (cf. Doi 
et al. 2005). 

3.2. Flux variability at A3 mm 

In NMA images at 100 GHz (A3 mm) of a field of 
view of ~ 1', which is roughly two times smaller than 
the JCMT image area, only a point-like single emission 
component was found at the nucleus. We detected the 
emission for 9 out of 15 epochs with NMA's sensitivities. 
Flux densities were measured in the image domain by 
elliptical Gaussian profile fitting with the JMFIT task of 
the MPS. We determined total errors in the flux measure- 
ments from the root sum square of errors of the Gaussian 
fitting and the flux scaling (15%). In our previous study 
for many LLAGNs (Doi ct al. 2005), we used the flux den- 
sity averaged over several of these epochs for NGG 4258. 

Even when the six undetected cases are allowed, the 
millimeter continuum emission tends to be significantly 
stronger than the centimeter emissions. The upper limit 
of dust contamination was estimated to be <2.2 mJy 



at 100 GHz, within a 15.1" region at the nucleus (§ 3.1). 
The upper limit of the contamination from low- frequency 
anomalous arms was estimated to be <0.5 mJy at 
100 GHz on the basis of a spectral extrapolation using a 
flux density of 3.4 mJy at 4.9 GHz for the nuclear region 
of 14" (Hummel et al. 1989) and a spectral index of —0.65 
(Hyman et al. 2001). Therefore, the contamination both 
from dust and anomalous arms was almost negligible for 
our NMA photometry with beams of ~3-7". 

Flux variability was detected at greater than 95% con- 
fldence. The variability showed no correlation with the 
array configuration (beam size), indicating that the de- 
tected variability did not result from a resolution effect. 
The NMA flux decreased from 18.1 ± 3.5 mJy at maxi- 
mum to 4.6 ± 1.7 mJy at minimum in three years, indi- 
cating L <0.7 (pc). A flux increase from 10.0 ± 1.3 mJy 
on December 21, 2001 to 18.1 ± 3.5 mJy on March 21, 
2002 (three months) results in L <0.13 (pc). We derive 

the size by the equation, L < c\ | ^ , where L is size, 
Sv is flux density, and t is observation time difference, as- 
suming an exponential light curve (Burbidge et al. 1974; 
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Fig. 2. — Contour map of the nuclear region in NGC 4258 
at 850 /im. Contour levels arc (7"'™^x (4, 5, 6, 7, 8, 9), where 
Q-rms — 5 mjy beam^^ is R.M.S. of image noise. Cross symbol 
represents the position of the NGC 4258 nucleus, and its positional 
accuracy is described by the length of the cross. Beam size, 15.1", 
is described in the lower left-hand corner of the image. 

Valtaoja et al. 1999). Slight diffuse contamination and 
the variability strongly demonstrate that almost all NMA 

fluxes originated in only a nuclear component rather than 
in any galactic components. Hence, our NMA photom- 
etry is substantially responsible for the active nucleus. 
This is the first insight into the nuclear radio emission at 
a millimeter band for NGC 4258. 

3.3. Spectral evolution at VLA wavelengths 

In VLA images at 5-43 GHz (A6-0.7 cm), a single 
emission component was identified at the nucleus in all 
cases we analyzed. Flux densities were measured with 
the JMFIT task of the MPS. The total errors were deter- 
mined from the root sum square of errors of the Gaussian 
fitting and the flux scaling. In all cases, emissions were 
almost unresolved: they were not extended beyond half 
of the FWHM of the synthesized beams. To estimate the 
contribution of diffuse components to our photometry, 
we also made images with uniform weighting, resulting 
in a slightly smaller beam (but in a less sensitive image). 
The resultant ratios of the measured flux densities with 
uniform weighting to that with natural weighting were 
between 0.80 and 1.10: 0.93 on average. The ratios indi- 
cate that almost all fluxes were from a region of < 0.1", 
corresponding to <3.6 pc at a distance of 7.2 Mpc. We, 
therefore, successfully carried out photometry without a 
significant resolution effect for the nucleus with the VLA- 
A and -B beams. 

The VLA spectrum for each epoch is shown in Fig. 4. 
Variability is apparently seen, compared to a power- 
law spectrum {a = 0.32) represented by dashed lines, 
which is determined by the least-square fit to all the data 
from 5 to 22 GHz. For each epoch, we also derived a 
power-law spectral index by least-square fitting to quasi- 
simultaneously obtained data from 5 to 22 GHz (col- 
umn [10] in Table 2), restflting in from a = 0.06 to 
a = 0.64; the average is (a) = 0.21 ± 0.04. The vari- 
ation of indices is significant at greater than 99.99% con- 



fidence. Note that the fitted power-law spectral indices 
remain inverted, a > 0, in all the epochs. 

The spectral features and their evolutions were appar- 
ently complicated. Several data points show significant 
departures from a simple linear trend. Highly inverted 
spectra are locally observed, for example a > 0.7 at 8.4- 
15 GHz on January 04, 1997 and September 05, 1998 
and at 15 22 GHz on January 05, 1997. On the other 
hand, fiat and marginally steep spectra are also locally 
observed. This indicates that the spectrum is made from 
multiple components. The first six epochs show spec- 
tral features rapidly changing in only a month, which 
suggests that the sizes of emitting sources shotild be re- 
stricted by these light-crossing times. There is a gradual 
increase at 22 GHz, from 2.2 ± 0.2 to 3.6 ± 0.4 in two 
weeks from December 22, 1996 to January 05, 1997, de- 
riving L <0.02 pc (<4000 AU). In addition, a fiux de- 
creases at 22 GHz in 13 days between March 06 and 19, 
1997, resulting in L < 0.02 pc. 

The trend of the inverted spectrum might continue into 
higher frequencies, which is inferred from the larger flux 
densities at 100 GHz at several epochs, although we have 
no simultaneous observation with VLA and NMA. A two- 
point spectrum between 15 GHz (A3. 6 cm) and 43 GHz 
(A7 mm), a = 0.59 ± 0.30, observed on December 21, 
2000 shows a sign that the trend of an inverted spectrum 
continues into millimeter wavelengths. The upper limit 
at 347 GHz (A850 /xm) from the JCMT data is too high 
to constrain whether the trend of the inverted spectrum 
continues into submillimeters. 

3.4. Comparisons to previous radio observations for 

NGC 4258 

3.4.1. Flux variability 

We discovered variable 100- GHz (A3mm) nuclear emis- 
sions at the nucleus. These millimeter emissions seem 
to be signiflcantly brighter than those at centimeters, 
which can also be seen in 43 GHz (A7mm) flux densities, 
although only two epoch data points are presented (Ta- 
ble 2). These millimeter emissions cannot represent the 
dust contribution because of the variability and the up- 
per limit of the expected dust spectrum based on the 
SCUBA/ JCMT image (§3.1). 

The significant flux variability at 22 GHz has been pre- 
viously suggested: the continuum component in VLBA 
images seems to vary several times (6'22GHz = 0.7- 
4.0 mJy) during 15 epochs in the period 1997-2001, al- 
though these measurements are based on a limited range 
of baseline lengths (Argon et al. 2007). Herrnstein et 
al. (1997) also showed flux variability of the continuum 
component at 22 GHz (from < 1 mJy to 6 mJy, based 
on both VLA and VLBA), which are correlated with the 
averaged flux density of maser emission in a systemic ve- 
locity range; the continuum emission is variable by up 
to 100% over timescales of a few wcciks. This range of 
flux densities and timescale of variability are consistent 
with our measurements at 22 GHz (§ 3.3). We have defi- 
nitely confirmed the presence of continuum flux variabil- 
ity in the NGC 4258 nucleus at high frequencies, even at 
100 GHz. 

3.4.2. Continuum spectra 
Ho & Ulvestad (2001) reported a slightly resolved core 
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Fig. 3. — Results of flux monitor for the NGC 4258 nucleus at 3 mm. The three plots on May 2005 arc from measurements at 95.725 GHz; 
the others are from those at 100.777 GHz. Symbols with arrows indicate 3cr upper limits. Flux variation is significant at greater than 95% 
confidence. 




Fig. 4. — Radio spectra for the NGC 4258 nucleus at 5-43 GHz obtained with the VLA A- or B-array. Dashed lines represent a fitted 
linear spectrum (5ighz = 1-2 mjy, a = 0.32, where Si, = S'ighz'^''^") for a^U data from 5 to 22 GHz. 
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of 5'i.4GHz ~ 2.73 and 55ghz = 2.18 (niJy), a ~ —0.18, in 
matched angular resolutions of ~ 1", in agreement with 
previous observations at a similar resolution (Vila et al. 
1990); Hyman et al. (2001) reported varying spectral in- 
dices between —0.18 and —0.41 between 1.4 and 5 GHz, 
i.e., slightly steep spectra, at the nucleus (on the other 
hand, the spectral index map shows relatively uniform 
distribution of a steep spectrum of a = —0.65 ± 0.10 
throughout "the anomalous arms"). An inverted spec- 
trum had been seen at higher frequencies: S^guz = 2.4 
and SisGHz = 3.2 (mJy), a = +0.4 in matched resolu- 
tions of - 3" (Turner '& Ho 1994). Thus, at 1.4 GHz, 
the extended emission of "the anomalous arms" possi- 
bly contaminates a nuclear flux density even with the 
A-configuration. Our observations at frequencies higher 
than 5 GHz of an extracted core emission (see, § 3.3) 
show inverted spectral indices in terms of power law fit- 
ted spectra in aU epochs (Col. (10) in Table 2) and an 
inverted spectral feature based on flux densities averaged 
through our observation epochs. 

Not only flux densities but also VLA spectral features 
are significantly variable in our measmements (Col. (10) 
in Table 2). Furthermore, the VLA spectral features can- 
not be represented with a simple power-law spectrum 
model (Fig. 4). Previously, Hyman et al. (2001) also 
showed variable spectral indices (slightly steep) between 
1.4 and 5 GHz toward the nucleus; Nagar et al. (2001) 
reported a folded spectral feature (2.1, 1.9, and 2.8 mJy 
at 5, 8.4, and 15 GHz, respectively) from VLA observa- 
tions at matched resolutions of 0.5". Thus, the contin- 
uum spectrum of the nucleus of NGC 4258 is definitely 
inverted, but not simple or stable. The physical origin of 
the spectral evolution is discussed in § 4. 

4. DISCUSSION 

4.1. Inverted spectra and a turnover frequency 

In this study, inverted radio spectra at centimeter 
wavelengths, probably up to millimeter wavelengths, 
were found for NGC 4258, and have properties similar 
to that found in Sgr A* and M 81. The fitted VLA 
spectra were always inverted but complicated and vari- 
able (§ 3.4.2), implying the presence of multiple emitting 
components. Such complicated and variable radio spec- 
tra are also seen in M 81 (Markofl" et al. 2008). On the 
other hand, the simultaneous spectrum at ~ 1-100 GHz 
of Sgr A* seems to be relatively simple with only a fold 
of the spectral index at ~ 10 GHz from a = 0.11 ± 0.03 
to a = 0.43 ± 0.04 (An et al. 2005); this is consistent 
with the result of other epochs (e.g., Falcke et al. 1998). 
We discuss later (§ 4.2) its physical origin as the outer 
region of nuclear components. 

The spectral turnover at a higher frequency, where the 
inverted spectrum becomes steep, cannot be determined 
for NGC 4258 using the available data, because of non- 
simultaneous observations among centimeter, millimeter, 
and submillimeter bands. The spectral turnover would 
provide important information about the innermost ra- 
dio emitting regions in the vicinity of the black hole. 
Turnover frequency is relevant to a radiiis where radia- 
tion becomes optically thin due to the extraction of ac- 
creting material via outflow (e.g., Di Matteo et al. 2000) 
or the inner edge of an accretion disk in the framework 
of accretion-flow models. In the framework of jet mod- 
els, the shock radius at the base of the jet is a major 



parameter that is responsible for a turnover frequency 
(Yuan et al. 2002). The observed spectral energy dis- 
tribution (SED) of Sgr A* suggests the presence of a 
turnover between the submillimeter and infrared regime, 
which has been explained by the radiation spectrum 
from the radius of the innermost edge of accretion flow 
(Narayan et al. 1995; Manmoto et al. 1997). The simul- 
taneous observations of M 81 at 100 GHz (A3mm) and 
230 GHz (Almm) present a turnover frequency of be- 
tween these frequencies (Schodel et al. 2007), which is 
in good agreement with the data presented by Renter & 
Lesch (1996); however the data were not acquired simul- 
taneously. The turnover frequency of M 81 is much lower 
than that of Sgr A*, suggesting that an effectively radi- 
ating site would be present at relatively outer region for 
M 81; the difference of Schwarzschild radii between them 
may be responsible for the difference of the turnover fre- 
quencies. On the other hand. Yuan et al. (2002) proposed 
that the SED of NGC 4258 can be modeled as a syn- 
chrotron jet dominated spectrum having a spectral peak 
in the infrared in order to satisfy both the weak radio 
and relatively strong infrared emissions (the infrared- 
optical spectral index indicates nonthermal emissions, 
a = -1.4±0.1; Chary et al. 2000) by the base of jets with 
a shock diameter of 5i?sch (-Rsch is the Schwarzschild ra- 
dius; Yuan et al. 2002). 

In the SED modeling for NGC 4258, observed radio 
emissions cannot be attributed to the same origin as the 
infrared peak by the base of jets; the outer component 
of nuclear jets is required (Yuan et al. 2002), because 
the emission of the base of jets that should come from 
the location of the dynamical center of the maser disk is 
suppressed to be < 220 ^Jy in VLBI images (Herrnstein 
et al. 1998). Thus, the determinations of turnover fre- 
quencies are crucial to reveal the physical origin of radio 
emissions from innermost regions. 

4.2. Physical origin of nuclear radio emissions 

The presence of jets has not been deflnitively confirmed 
or ruled out for Sgr A*. M 81 shows structures consisting 
of a core and jet-like elongation at 1-43 GHz (Bieten- 
holz et al. 2004; Ros & Perez-Torres 2008), and the core 
contributions dominate total radio flux in VLBI images. 
The location of the core (brightness peak) moves south- 
west with increasing frequency (1.7-15 GHz), and this 
trend continues up to the stationary position of a puta- 
tive black hole (Bietenholz et al. 2004). The line of sight 
toward the nucleus is opaque at these radio frequencies. 
Hence, the core emissions cannot be responsible for an 
accretion disk except due to jets. Thus, the observed in- 
verted spectrum in total flux definitively comes from jets 
for M 81 at these frequencies. The size of the emitting 
region inferred from the time scales of observed millime- 
ter variability is < 25i?sch (Schodel et al. 2007, see also 
Sakamoto et al. 2001), which is possibly consistent with 
the radio emitting size of the accretion disk predicted by 
the ADAF model. The spectral tmnovcr at 100 -230 GHz 
(Schodel et al. 2007) implies that millimetric observations 
could see through to the vicinity of the black hole. 

The NGC 4258 nucleus has a similar situation in terms 
of the physical origin of nuclear radio emission. Jet struc- 
tures whose brightness peaks are significantly shifted 
from a dynamical center are clearly seen in VLBI im- 
ages at 1.5 and 22 GHz (Cecil et al. 2000; Herrnstein et 
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al. 1997; Argon et al. 2007). The observed radio luminos- 
ity of NGC 4258 at 5 GHz is higher than the maximal 
luminosity expected from the ADAF model and black 
hole mass, implying that the radio emission is predomi- 
nantly from jets, although its optical luminosity could be 
provided by the ADAF (Wu & Cao 2005); M 81 is also 
in the same situation. Hence, the observed VLA spec- 
tra and their time evolution in our study for NGC 4258 
should be caused by jets; we discuss the spectral evolu- 
tion in detail in § 4.3. 

The emission from the location of the dynamical cen- 
ter is < 220 AiJy at 22 GHz (Herrnstcin ct al. 1998), 
implying higher frequencies are needed to see through. 
A spectral turnover at a higher frequency, where the 
source is expected to become optically thin, cannot be 
determined using available data. The observed millime- 
ter variability during three months in our study suggests 
a size of < 0.17 pc (§ 3.2), which is too large to give a 
constraint for the emission from the scale of the accre- 
tion disk. Multi-frequency millimeter observations aim- 
ing to investigate intra-day spectral evolution are crucial 
to determine a spectral turnover and the contribution of 
emission from the scale of the accretion disk. 

4.3. An interpretation of observed spectral evolution 

Radio light curves of Sgr A* at 22 and 43 GHz for 
a low level of flare activity can be interpreted by con- 
sidering the ejection and adiabatic expansion of a self- 
absorbed plasma blob, showing a spectral evolution that 
a spectral peak frequency shifts toward lower frequencies 
(Yusef-Zadeh et al. 2006). In this scheme, spectral evo- 
lution involving multiple blobs is predicted to propagate 
from a high frequency to lower frequencies. Similar ap- 
parent spectral changes are also seen in M 81 (Markoff et 
al. 2008). On the contrary, the observed VLA spectra of 
NGC 4258 seem to include an apparent spectral evolution 
of a flare-up propagating from 8.4 to 22 GHz in the first 
four epochs (Figure 4) , which cannot be explained by the 
simple model of an expanding blob; this event began in 
the period December 07-22, 1996 and ended in the pe- 
riod December 29-31, 1996, i.e., during 7 <t <24 days. 
In this paper, for NGC 4258, we propose the model of a 
series of self-absorbed blobs affected by external absorp- 
tion (free-free absorption), as described below. 

The intensity peak of jets of NGC 4258 at 22 GHz is 
located about 0.4 mas (corresponding to 0.014 pc) north 
of the dynamical center of the disk, which is viewed al- 
most edge-on (inclination of 82°; Herrnstcin et al. 1999). 
The relative faintness of the southern jet component may 
arise because of free-free absorption due to an interven- 
ing layer of ionized atomic gas above the molecular disk 
(Herrnstcin et al. 1997). The emission from the loca- 
tion of the dynamical center is suppressed to < 220 /iJy, 
which has been discussed as an ADAF-emitting radius 
smaller than lOOiigch (Herrnstcin et al. 1998) or as the 
opaque nozzle of a relativistic plasma flow (Blandford & 
Konigl 1979; Herrnstcin et al. 1997). In this paper, we 
consider another case that possibly the disk inflates to be 
thickc;r due to higher temperatures at the inside of the 
maser orbits and could be an absorber toward the root 
of the northern jet in addition to the nucleus and the 
root of the southern jet. The thickness of outer (0.13- 
0.26 pc in radii) maser disk is 5 /xarcsec (40 AU) (Argon 
et al. 2007); we consider the thickness of the putative 



inner (< 0.13 pc) disk to be - 0.2 mas 1400 AU). 
We assume that the observed spectral evolution is a 
composite of a series of emerging synchrotron emitting 
blobs of spectral-peak flux densities of ~ 0.4, ~ 0.5, and 
~ 0.7 mJy at 8.4, 15 and 22 GHz, respectively; these 
components are supposed to provide residual fluxes be- 
tween December 07-31, 1996. Under the equipartition 
condition between the energy densities of particles and 
magnetic fields in a uniform spherical plasma, the diam- 
eters of these blobs are - 90, ~ 50, and ~ 40 AU (~ 110, 
~ 70, and ~ 60i?sch) at 8.4, 15 and 22 GHz, respectively; 
the total extent of the jet consisting of three blobs would 
be > 180 AU unless they overlapped each other. If these 
blobs get away from the absorption region while advanc- 
ing north from the dynamical center, in order of the 8.4-, 
15-, and 22-GHz blobs, then the observed spectral evo- 
lution propagating from 5 to 22 GHz can be explained. 
The advancing speed should be > 0.044c (c is the speed 
of light), because of the total extent of jets of > 180 AU 
and a duration of the event of < 24 days. The flare up 
of 22 GHz seems to begin before December 22 and to 
end after December 29; the duration of 22-GHz flaring 
up seems to be > 7 days, which means an advancing 
speed of < 0.12c to pass over the diameter of the 22- 
GHz blob. Consequently, the putative jet speed would 
be sub-relativistic, 0.044c-0.12c. 

Our model requires no spectral evolution of each in- 
dividual blob during the event, i.e., no adiabatic expan- 
sion. Using the model of an adiabatically expanding syn- 
chrotron blob (van der Laan 1966) to avoid the decrease 
of spectral peak frequency from 15 to 8.4 GHz during the 
event of 7 < t < 24 days, an expansion speed should be 
limited to 0.002c < v < 0.006c, which is much smaller 
than the putative bulk speed. Thus, our model described 
above is valid only in a dynamical situation dominated 
by the time scale of bulk advancing jets rather than by 
expansion. 

A moving absorber in front of emitting blobs is an- 
other idea to explain the spectral evolution propagating 
to higher frequencies. A Kepler rotation velocity of the 
maser disk is 1140 km s~^ (0.0038c) at the radius of the 
innermost maser disk of 0.13 pc (e.g., Miyoshi et al. 1995; 
Argon et al. 2007). To provide the rotation speed of the 
required > 0.044c, a rotation radius should be < 200 AU, 
which is close to the same as or less than the size of the 
emitting region. Thus, a moving absorber cannot be re- 
sponsible for the observed spectral evolution. 

4.4. NGC 4258 another scaled-up version of Sgr A * 

We found inverted radio spectra in the NGC 4258 nu- 
cleus, analogous to Sgr A* and M 81. We discuss the 
analogy and differences on the NGC 4258 as another 
scaled-up version of Sgr A*, in addition to M 81, in ref- 
erence to Table 1. 

The NGC 4258 is a radio source with a linear scale 
800 times larger (80 times larger in i?sch unit) than that 
of Sgr A*. The black hole mass of NGC 4258 is up- 
scaled with reference to Sgr A*, in the same order as that 
of M 81. However, several differences are seen between 
NGC 4258 and M 81. 

The X-ray luminosity of NGC 4258 is flve times larger 
than that of M 81, and its 10 times higher Lx/^Edd, 
which is thought to be proportional to the Eddington ra- 
tio iboi/^Edd that is deflned as the ratio of bolometric 
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to Eddington luminosities, indicates a higher accreting 
system compared to M 81. However, the radio power is 
10 times smaller. According to "a fundamental plane of 
black hole activity" (Merloni et al. 2003, see also Falcke 
et al. 2004), radio luminosity seems to depend on the re- 
lationship of LsGHz = lO'^-^^L'^'^M^-l^, where the 5-GHz 
core luminosity Lsghz and the 2-10 kcV X-ray luminos- 
ity Lx in ergs s~^, and the black hole mass Mbh in Mq, 
with a scatter of cr(I/5GHz) = 10°-^®. The radio lumi- 
nosities the relation predicts are ~ 1/100 and 1/5 of the 
observed ones for NGC 4258 and M 81, respectively, as- 
suming a ~ 0.2 at 5-22 GHz. Therefore, the radio emis- 
sion of the NGC 4258 nucleus is very weak. On the other 
hand, the linear scale of the radio source is four times 
larger than that of M 81 (also eight times larger scaled 
in i?sch unit). Hence, NGC 4258 has significantly lower- 
brightness jets. The core shift of NGC 4258 is eight times 
larger than the expected core shift of ~ 0.09 mas at 
22 GHz for M 81^. This is contrary to commonsense: 
a core shift tends to be larger in a higher radio luminos- 
ity. The difference of inclinations {i = 82° for NGC 4258; 
Miyoshi et al. 1995; Herrnstein et al. 1999, and i = 14° 
for M 81; Devereux et al. 2003) can account for half of 
the factor of 8, and thus another parameter providing a 
large factor may be required. A possible parameter is 
the jet Lorentz factor 7 (or jet speed /3); ^2.5-3 times 
smaller /3 could be responsible for the larger core shift 
of NGC 4258, considering the radio luminosities and in- 
clinations under a sub-relativistic regime (Eq. (11) in 
Lobanov 1998). Thus, compared to that of M 81, the 
jets of NGC 4258 may be less accelerated, which is con- 
sistent with the pictiue of significantly lower-brightness 
jets as pointed out above. Otherwise, free-free absorp- 
tion would contribute the larger core shift of NGC 4258. 
Such a situation can be relevant to the discussion of the 
VLA spectral evolution with an absorber at the inner 
disk (§ 4.3). 

We discussed above the differences of the resolved jets 
of NGC 4258 and M 81. On the other hand, the presence 
of jets has not been definitively confirmed or ruled out for 
Sgr A*, and the physical scales of radio emitters are quite 
different (a factor of ^ 1000 in linear scale, or ~ 100 in 
Rsch)- Although Sgr A* is a ^ 10^ times lower Edding- 
ton ratio system compared to NGC 4258 and M 81, it is 
in good agreement with the fundamental plane (the ex- 
pected radio luminosity is 6.4 x 10^^ W Hz~^ at 5 GHz). 
This is suggestive of the jet origin in a radiatively inef- 
ficient accreting system (Merloni et al. 2003). However, 
Sgr A* seems to have no ability to emanate jets up to 
kpc scales, unlike NGC 4258 (anomalous arms; Hyman 
et al. 2001) and M 81 (Kaufman et al. 1996). Such great 
distinctions of physical scales of jets might be relevant to 
the difference between the relatively simple radio spec- 
tra of Sgr A* and complicated spectra of NGC 4258 and 
M 81. Thus, in addition to M 81, NGC 4258 is possibly 
an up-scaled version of Sgr A* and shows some appar- 
ently difii'erent radio properties, which might be related 
to physical scale. 

5. SUMMARY 



We investigated the radio properties of the contin- 
uum emission from the nucleus of the nearby LLAGN 
NGC 4258, and discussed analogies to Sgr A* and M 81 
as a scaled- up version of Sgr A*. The summaries are 
described below: 

• We discovered a large amplitude variable emission 
at 100 GHz (A3mm) using the NMA. The emission 
was detected at 9 out of 15 epochs. The detected 

flux densities were higher than those at centime- 
ters (not simultaneous). Signiflcant contamination 
from dust emission is ruled out by the variability 
and the estimation using JCMT data at 345 GHz. 

• Quasi-simultaneous multi-frequency observations 
using the VLA at 5-22 GHz showed inverted spec- 
tra in all ten epochs. The spectra were compli- 
cated and variable, indicating the presence of mul- 
tiple blobs in nuclear jets. The observed spectral 
evohition cannot be explained by the model of an 
adiabatically expanding blob that is expected to 
show a spectral peak propagating to lower frequen- 
cies, while a propagation to higher frequencies was 
apparently seen in less than one month. We pro- 
posed partially absorbed advancing components to 
explain the spectral evolution. 

• The observed radio property of a highly variable 
and inverted spectrum for the NGC 4258 nucleus 
is similar to that of Sgr A* and M 81, even with 
more than three orders of magnitude larger radio 
power, 10 times larger black hole mass, and a more 
than six orders of magnitude higher Eddington ra- 
tio compared to Sgr A. In addition to the LLAGN 
M 81, NGC 4258 is possibly another up-scaled ver- 
sion of Sgr A*. 
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^ On the basis of astrometric measurements at 1.7-15 GHz, the etenholz et al 2004) 

stationary position of a putative black hole coincides with the 50% 
of contour levels of brightness distribution of each frequency (Bi- 



REFERENCES 



An, T., Goss, W. M., Zhao, J.-H., Hong, X. Y., Roy, S., Rao, A. R, 

& Shcn, Z.-Q. 2005, ApJ, 634, L49 
Argon, A. L., Greenhill, L. J., Rcid, M. J., Moran, J. M., & 

Humphreys, E. M. L. 2007, ApJ, 659, 1040 
Baganoff, F. K., ct al. 200.3, ApJ, 591, 891 
Bietenholz, M. F., ct al. 1996, ApJ, 457, 604 

Bietenholz, M. F., Bartcl, N., & Rupcn, M. P. 2004, ApJ, 615, 173 
Blandford, R. D., & Konigl, A. 1979, ApJ, 232, 34 
Bower, G. C., Wright, M. C. H., Falcke, H., & Backer, D. C. 2003, 
ApJ, 588, 331 

Bower, G. C., Goss, W. M., Falcke, H., Backer, D. C., & Lithwick, 

Y. 2006, ApJ, 648, L127 
Burbidge, G. R., Jones, T. W., & Odell, S. L. 1974, ApJ, 193, 43 
Cecil, G., et al. 2000, ApJ, 536, 675 

Chary, R., Becklin, E. E., Evans, A. S., Neugebauer, G., Scoville, 
N. Z., Matthews, K., & Ressler, M. E. 2000, ApJ, 531, 756 

Devereux, N., Ford, H., Tsvetanov, Z., & Jacoby, G. 2003, AJ, 125, 
1226 

Di Mattco, T., Quataert, E., Allen, S. W., Narayan, R., & Fabian, 

A. C. 2000, MNRAS, 311, 507 
Doi, A., Kameno, S., Kolmo, K., Nakanislii, K., & Inoue, M. 2005, 

MNRAS, 363, 692 
Dunne, L., Eales, S., Edmunds, M., Ivison, R., Alexander, P., & 

Clements, D. L. 2000, MNRAS, 315, 115 
Eisenhauer, F., Schodel, R., Genzel, R., Ott, T., Tecza, M., Abuter, 

R., Eckart, A., & Alexander, T. 2003, ApJ, 597, L121 
Falcke, H., Goss, W. M., Matsuo, H., Tcuben, P., Zhao, J.-H., & 

Zylka, R. 1998, ApJ, 499, 731 
Falcke, H., Kording, E., & Markoff, S. 2004, A&A, 414, 895 
Falcke, H., & Bicrmann, P. L. 1999, A&A, 342, 49 
Frcedman, W. L., ct al. 2001, ApJ, 553, 47 

Fruscionc, A., Greenhill, L. J., Filippenko, A. V., Moran, J. M., 

Herrnstcin, J. R., & Gallc, E. 2005, ApJ, 624, 103 
Gillesscn, S., Eisenhauer, F., Trippe, S., Alexander, T., Genzel, R., 

Martins, F., & Ott, T. 2009, ApJ, 692, 1075 
Heifer, T. T., Thornley, M. D., Regan, M. W., Wong, T., Sheth, 

K., Vogel, S. N., Blitz, L., & Bock, D. C.-J. 2003, ApJS, 145, 259 
Herrnstein, J. R., Moran, J. M., Greenhill, L. J., Diamond, P. J., 

Miyoshi, M., Nakai, N., & Inoue, M. 1997, ApJ, 475, L17 
Herrnstein, J. R., Greenhill, L. J., Moran, J. M., Diamond, P. J., 

Inoue, M., Nakai, N., & Miyoshi, M. 1998, ApJ, 497, L69 
Herrnstein, J. R., et al. 1999. Nature, 400, 539 
Hildebrand, R. H. 1983, QJRAS, 24, 267 

Ho, L. C, Filippenko, A. V., & Sargent, W. L. W. 1997, ApJS, 
112, 315 

Ho, L. C, & Ulvestad, J. S. 2001, ApJS, 133, 77 
Holland, W. S., et al. 1999, MNRAS, 303, 659 
Hummel, E., Krause, M., & Lesch, H. 1989, A&A, 211, 266 
Hyman, S. D., Calle, D., Weiler, K. W., Lacey, C. K., Van Dyk, 
S. D., & Sramek, R. 2001, ApJ, 551, 702 



Ishisaki, Y., ct al. 1996, PASJ, 48, 237 

Kaufman, M., Bash, F. N., Crane, P. C, & Jacoby, G. H. 1996, 

AJ, 112, 1021 
Lobanov, A. P. 1998, A&A, 330, 79 

Manmoto, T., Mineshige, S., & Kusunose, M. 1997, ApJ, 489, 791 
Mao, R.-Q., Schulz, A., Henkel, C, Mauersberger, R., Muders, D., 

& Dinh-V-Trung 2010, ApJ, 724, 1336 
Markoff, S., et al. 2008, ApJ, 681, 905 

Merloni, A., Heinz, S., & di Matteo, T. 2003, MNRAS, 345, 1057 
Miyoshi, M., Moran, J., Herrnstein, J., Greenhill, L., Nakai, N., 

Diamond, P., & Inoue, M. 1995, Nature, 373, 127 
Nagar, N. M., Wilson, A. S., & Falcke, H. 2001, ApJ, 559, L87 
Nakai, N., Inoue, M., & Miyoshi, M. 1993, Nature, 361, 45 
Narayan, R., Yi, I., & Mahadevan, R. 1995, Nature, 374, 623 
Ozel, F., Psaltis, D., & Narayan, R. 2000, ApJ, 541, 234 
Regan, M. W., Thornley, M. D., Heifer, T. T., Sheth, K., Wong, 

T., Vogel, S. N., Blitz, L., & Bock, D. C.-J. 2001, ApJ, 561, 218 
Reid, M. J. 1993, ARA&A, 31, 345 
Renter, H.-R, & Lesch, H. 1996, A&A, 310, L5 
Ros, E., & Perez-Torres, M. 2008, The role of VLBI in the Golden 

Age for Radio Astronomy 
Sakamoto, K., Fukuda, H., Wada, K., & Habe, A. 2001, AJ, 122, 

1319 

Sawada-Satoh, S., Ho, P. T. P., MuUer, S., Matsushita, S., & Lim, 

J. 2007, ApJ, 658, 851 
Shen, Z.-Q., Lo, K. Y., Liang, M.-C, Ho, P. T. P., & Zhao, J.-H. 

2005, Nature, 438, 62 
Schodel, R., Krips, M., Markoff, S., Neri, R., & Eckart, A. 2007, 

A&A, 463, 551 

Tsutsumi, T., Morita, K.-L, & Umeyama, S. 1997, Astronomical 
Data Analysis Software and Systems VI, 125, 50 

Turner, J. L., & Ho, P. T. P. 1994, ApJ, 421, 122 

Valtaoja, E., Lahteenmaki, A., Tcrasranta, H., & Laincla, M. 1999, 
ApJS, 120, 95 

van der Kruit, P. C, Oort, J. H., & Mathcwson, D. S. 1972, A&A, 
21, 169 

van der Laan, H. 1966, Nature, 211, 1131 

Vila, M. B., Pedlar, A., Davies, R. D., Hummel, E., & Axon, D. J. 

1990, MNRAS, 242, 379 
Wu, Q., & Cao, X. 2005, ApJ, 621, 130 

Yamada, S., Itoh, T., Makisliima, K., & Nakazawa, K. 2009, PASJ, 
61, .309 

Yuan, F., Markoff, S., & Falcke, H. 2002, A&A, 383, 854 
Yuan, F., Quataert, E., & Narayan, R. 2003, ApJ, 598, 301 
Yusef-Zadeh, F., Roberts, D., Wardle, M., Heinke, C. O., & Bower, 
G. C. 2006, ApJ, 650, 189 



